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We report measurements of branching fractions and time-dependent CP asymmetries in B° — > 
D + D~ and B° — > D* ± D^ 1 decays using a data sample that contains (772 ± 11) x 10 6 BB pairs 
collected at the T(4S) resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. We determine the branching fractions to be B (B° -> D + D~) = (2.12 ± 0.16 ± 0.18) x 10" 
and B (B° -> D f± D T ) = (6.14 ± 0.29 ± 0.50) x 10~ 4 . We measure CP asymmetry parameters 
S D+D - = -1.06i^4±0.08andC D+D - = -0.43±0.16±0.05 in B° -> D + D' &m\A D -D = +0.06± 
0.05±0.02, Sd-d = -0.78±0.15±0.05, C d *d = -0.01±0.11±0.04, AS d -d = -0.13±0.15±0.04 
and ACd*d = +0.12 ± 0.11 ± 0.03 in B° — > D* ± D^, where the first uncertainty is statistical and 
the second is systematic. We exclude the conservation of CP symmetry in both decays at equal to 
or greater than 4a significance. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Ff 
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In the standard model (SM) of electroweak interac- 
tions, the effect of CP violation is explained by a single 
complex phase in the three-family Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix p]. Both the 
Belle and BaBar Collaborations experimentally estab- 
lished this effect [21 [3] and precisely determined the 
parameter sin2^i by measurements of mixing-induced 
CP asymmetries in b — > (cc)s transitions, where 4>\ — 

™&l-v cd v; b /v td v t * b ] m- 

In b — > ccd transitions such as B° — > D^^D^ decays, 
the dominant contributions are Cabibbo-disfavored but 
color-allowed tree-level diagrams and the corresponding 
mixing-induced CP asymmetries are directly related to 
sin 20i. In addition, b — > d penguin diagrams that may 
have different weak phases can contribute to these de- 
cays. Theoretical considerations based on models using 
factorization approximations and heavy quark symmetry 
predict the corrections to mixing-induced CP violation 
to be a few percent and possible direct CP violation to 
be negligibly small [7j. 

CP violation in b — > ccd transitions has been studied 
previously by the Belle and BaBar Collaborations. In 



B° -» D + D~ decays using a data sample of 535 x 10 6 BB 
pairs, Belle found evidence of a large direct CP viola- 
tion: C D + D - = —0.91 ± 0.23 ± 0.06 corresponding to a 
3.2(7 deviation from zero [8j[9], in contradiction to the- 
oretical expectations [7]. This deviation was not con- 
firmed by BaBar and has not been observed in other 
B° D*±Z?(*)=F decay modes [Mill- 
In this article we present measurements of branching 
fractions and CP violating asymmetries in the decays 
B" ->■ D+D- and B° -> D^D^ using the final data 
sample of the Belle experiment. 

The decay rate of a neutral B meson decaying to a CP 
eigenstate such as D + D~ is given by 

e -|At|/T B o 

f D+D - (At) = — {1 + q[S D+D - sm(Am d At) 

4r B o 

-C D+D -cos(Am d At)}}, (1) 

where q = +1 (—1) represents the 6-flavor charge when 
the accompanying B meson is tagged as a B° (B°), and 
At represents the proper time interval between the two 
neutral B decays in an T (AS) event. The B° lifetime 
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is denoted by Tgo and the mass difference between the 
two neutral B mass eigenstates by Am^. The parameters 
Sd+jj- and Cjj+jj- measure mixing-induced and direct 
CP violation, respectively [5]. 

Unlike D + D~ , D*+ D~ and D*~ D+ are not CP eigen- 
states. The decay rate of neutral B mesons decaying to 
these states has four flavor-charge configurations and can 
be expressed as [Ml IT4] 

-|A*|/t b o 

f D ,± D *(At) = (1 ± = 

8r B o 

x{l + q[(S D . D ±AS D . D ) sm(Am d At) 

-{C D , D ±AC D . D )cos(Am d At)]}, (2) 

where the + (— ) sign represents the D* + D~ (D*~D + ) 
final state. The time- and flavor-integrated charge asym- 
metry Ad* d measures direct CP violation. The quantity 
3d* d parameterizes mixing-induced CP violation and 
Cd*d parameterizes flavor-dependent direct CP viola- 
tion. The quantities ACd*d and ASd*d are not sen- 
sitive to CP violation. The parameter ACd* d describes 
the asymmetry between the rates T(B° —> D*~D + ) + 
T(B° -> D*+D-) and T{B° -> D*+D-) + T(B° -> 
D*~D + ). The parameter ASd*d is related to the rel- 
ative strong phase between the amplitudes contributing 
to the decays. 

This analysis is based on a data sample containing 
(772 ± 11) x 10 6 BB pairs collected at the T(45) reso- 
nance with the Belle detector at the KEKB asymmetric- 
energy e + e~ collider [T5]. The T(45) is produced with 
a Lorentz boost of (3j — 0.425 close to an axis along the 
e~ beam, which allows the determination of At from the 
displacement of decay vertices of both B mesons. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that is described in detail in Ref. 16J. The 
present analysis uses for track reconstruction and parti- 
cle identification a silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like ar- 
rangement of time-of-flight scintillation counters (TOF) , 
and an electromagnetic calorimeter composed of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. 

Reconstructed charged tracks are required to have a 
transverse (longitudinal) distance of closest approach to 
the interaction point (IP) of less than 2 (4) cm. For iden- 
tification of charged particles (PID), measurements of 
specific energy loss in the CDC and measurements from 
the ACC and TOF are combined in an likelihood-ratio 
approach. The selection requirement on the combined 
PID quantity has a kaon (pion) identification efficiency 
of 91% (99%) with an associated pion (kaon) misiden- 
tification rate of 2% (18%). Charged tracks are also re- 
quired to be not positively identified as electrons by mea- 
surements of shower shapes and energy deposited in the 
ECL. Neutral pions are reconstructed from two photons 
detected in the ECL with each photon having an energy 
greater than 30 MeV. The invariant mass of the photon 



pair is required to be within 15 MeV/c 2 of the nomi- 
nal 7T° mass (corresponding to a width of 3.3ct). For w° 
candidates a kinematic fit to the IP profile with a mass 
constraint is performed. Neutral kaons are reconstructed 
in the decay mode Kg — > tt + it~. The invariant mass 
of the 7r + 7r~ pair is required to be within 15 MeV/c 2 
of the nominal K^ mass (5.8ct). Additional momentum- 
dependent selection requirements consider the possible 
displacement of the Kg decay vertices from the IP [T7j . 

Charged D mesons are reconstructed in the decay 
modes D+ -> K-j: + tt + and D+ -> K s n+ [18 . The 
invariant mass of D + candidates is required to be within 
12 MeV/c 2 of the nominal mass (3.4cr in D+ -> K~n + ir + 
and 2.9a in D + — > Kgir + ). Neutral D mesons are re- 
constructed in the decay modes D° — > K~n + , D° — > 
K-ir+ir+n-, K^it+tt- and D° -> K-jt+tt . The in- 
variant mass of D° candidates is required to be within 
15 MeV/c 2 (3.3<7 — 3.7a) of the nominal mass, except 
for the D° — > K~it + ir° decay mode where a require- 
ment of 32 MeV/c 2 (3. Oct) is applied. We reconstruct 
D* + mesons in the decay modes D* + — » D°ir + and 
D* + — > D + ir°. The momentum resolution of charged 
low momentum pions from D* + decays, referred to as 
soft pions, is improved by a kinematic fit in which the 
soft pion is constrained to the D* + decay vertex deter- 
mined from a kinematic fit of D candidates constrained to 
originate from the IP profile. The difference of invariant 
masses of the D* + and D° (D + ) candidates is required 
to be within 1.5 MeV/c 2 (2.5 MeV/c 2 ) (3.1a - 3.7ct) of 
the nominal mass difference, except for modes involving 
D° — > K~ir + ir + ir~ and D Q — > K~tt + it decays where a 
requirement of 2 MeV/c 2 (4.6ct and 3.2ct) is applied. 

Neutral B mesons are reconstructed by com- 
bining £)(*)+ and D~ candidates, and se- 
lected by the bcam-cncrgy-constrained mass 
M bc = v^boam/c 2 ) 2 - (Ps/c) 2 and the energy 
difference AE = E* B - E* cam , where E* cam is 
the energy of the beam and and Eg are the 
momentum and energy of the B° candidates in 
the center-of-mass frame (cm.). The selected re- 
gions are 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and 
-50 MeV < AE < 100 MeV. The lower boundary in 
AE was chosen to exclude reflections from misidentified 
B° -> Z>+I)(*5- decays that populate the M hc signal 
region at AE w -75 MeV. 

In B° -> D+D- (B° -> D* ± L> = F) ! a f ter applying the 
above selection requirements, 12% (16%) of the signal 
events contain more than one B° candidate. In this 
case the candidate with the smallest quadratic sum of 
deviations of reconstructed invariant masses of D daugh- 
ters (and mass differences of D* + daughters) from nom- 
inal values, divided by the width of corresponding signal 
peaks, is selected. This requirement selects the correct 
candidate with a probability of 96% (92%). 

In B° -> D+D- unlike in B° -> D^D^ the 
major source of background arises from e + e~ — > qq 
(q G {u,d, s,c}) continuum events. This background is 
suppressed by a neural network (NN) implemented by 
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the NeuroBayes package [TH] that combines information 
about the event topology. Observables included in the 
NN are cos 9* B , where 9* B is the polar angle of the B 
candidate with respect to the beam direction in the cm. 
frame, a combination of 16 modified Fox- Wolfram mo- 
ments |20j . and the momentum flow in nine concentric 
cones around the thrust axis of the B° candidate [2T] . 
The requirement on the NN selection rejects 64% of the 
background while retaining 92% of the signal. 

The signal yields are obtained by two-dimensional un- 
binned extended maximum likelihood fits to the Afb c and 
AE distributions. The Mb c distributions are parameter- 
ized by a Gaussian function for the signal component 
and by an empirically determined threshold function in- 
troduced by the ARGUS Collaboration [55] for the back- 
ground component. The AE distributions are parame- 
terized by the sum of two Gaussian functions (the sum of 
a Gaussian function and an empirically determined func- 
tion introduced by the Crystal Ball Collaboration [2"3"] ) 
with common mean for the signal component in B° — > 
D + D~ (B a — > D*+D+) and by a linear function for the 
background component. The shape parameters of signal 
components in B° D+D- (B° ->■ D*+D+) are fixed to 
values obtained from B° ->■ D+D- (B° D+D*-) data 
distributions, where the relative widths and fractions of 
the signal components in AE are fixed to values obtained 
from Monte Carlo (MC) simulation studies. The Mb c and 
AE distributions and fit projections are shown in Fig. [I] 
For B° D+D~ the obtained yields are 221.4±18.6 sig- 



nal events in the (K ir+n+)(K+Tr tt 



48.0 ± 8.9 signal events in the (K tt+^ 



final state and 
+)(K%tt-) final 



state. 

For B° D*+D+, we obtain a yield of 886.8 ± 39.3 
signal events in all reconstructed modes combined. Of 
these, the yield in modes involving D*+ — > D°tt+ decays 
only is 769.2 ± 36.0 signal events. 

Decays such as B° D^~K* + , B° D^-K tt+ 
and B° — > D^- n+ tt+ it- have the same final states 
as the reconstructed B° —> D+D^~ decay modes and 
can possibly populate the M\, c and AE signal region. 
The contributions of such decays, referred to as peak- 
ing background, are estimated from D mass sidebands 
and subtracted in the signal yields given above. For 
B° -> D+D- (B° -> D*+D+), we find a contribution of 
0.7 ± 1.5 (4.7 ± 2.1) peaking background events from fits 



to D — > K9.tt mass sidebands. The D 



K+TT~ 



mass sidebands are considered to be free of peaking back- 
ground and no background subtraction is performed. 
This assumption has been tested by MC simulations and 
no peaking background is found in the data sidebands. 

The reconstruction efficiencies are obtained from MC 
simulations of signal decays and have been corrected to 
account for PID selection efficiency differences between 
MC simulations and data. To exclude systematic ef- 
fects in the determination of reconstruction efficiencies 
associated with soft neutral pions, only modes involving 
D*+ -> D°n+ decays are used in the B° -> D*+D+ 
branching fraction measurement. 




5.20 5.22 5.24 5.26 5.28 5.30 

(C) M bc (GeV/e 2 ) 



FIG. 1: Afbc and AE distributions (data points with error 
bars) and fit projections (solid lines) for (a)-(b) B° —¥ D+ D~ 
and (c)-(d) B° -> D*+D+ decays . The dotted (dashed) lines 
represent projections of signal (background) fit components. 
A \AE\ < 30 MeV (M bc > 5.27 GeV/c 2 ) requirement is ap- 
plied in plotting the Mb c (AE) distributions. 



The branching fractions are calculated from signal 
yields, reconstruction efficiencies, the number of BB 
events and current world averages of D°, D+ and D*+ 
branching fractions [21] • The branching fraction for 
B° — » D+D- decays is calculated as the weighted aver- 
age of the branching fractions determined for each of both 
reconstructed decay modes separately. The branching 
fraction for B° — > D*+D+ decays is determined by the 
signal yield in all modes and the average reconstruction 
efficiency weighted by the D branching fractions. The 
determined branching fractions are B [B° — s- D+D~) = 
(2.12 ±0.16 ±0.18) x 10" 4 and B(B°^D* ± D+) = 
(6.14 ± 0.29 ± 0.50) x 10~ 4 . The systematic uncertainties 



TABLE I: Summary of systematic uncertainties of the B° 
D+D- and B° -> D*+D+ branching fractions (in %). 



Source 


D+D- 


D* ± D+ 


Track reconstruction efficiency 


2.0 


4.1 


Kg reconstruction efficiency 


0.7 


0.7 


7T° reconstruction efficiency 




1.6 


K/n selection efficiency 


5.5 


5.3 


Event reconstruction efficiency 


1.0 


0.1 


Continuum suppression 


4.1 




Fit models 


1.1 


0.6 


D branching fractions 


4.3 


3.9 


Number of BB events 


1.4 


1.4 


Total 


8.6 


8.1 



5 



of the measured branching fractions are summarized in 
Table|l] The uncertainties due to track, Kg and 7r° recon- 
struction efficiency and the uncertainty due to the K/tt 
selection efficiency have been estimated using studies of 
D*+ decays with MC simulations and data. The effect 
on the event reconstruction efficiencies due to broader D 
mass distributions for data and the corresponding selec- 
tion is studied by a MC/data comparison and assigned as 
a systematic uncertainty. As the systematic uncertainty 
of the applied continuum suppression in B° —> D+D~, 
the maximum variation of signal yields in a MC/data 
comparison of the neural networks using B° — > D+D~ 
decays is assigned. The contributions due to the fit mod- 
els are estimated by varying the fixed parameters within 
their uncertainties. The contributions due to uncertain- 
ties of the D°, D+ and D*+ branching fractions and of 
the number of BB events are obtained by propagation 
of the appropriate uncertainties. The total systematic 
uncertainties are obtained by adding all contributions in 
quadrature. 



1 



1 




FIG. 2: Top: At distributions (data points with error bars) 
of (a) B° -> D+D~ and (b) B° -> D*+ D~ + B° -> D*~ D+ 
candidates associated with high quality flavor tags (r > 0.5). 
The lines show projections of the sum of signal and back- 
ground components in the fit. The signal purity for r > 0.5 
is 69% (66%) for B° -> D+D' {B° -> D* ± D+). Bottom: 
The CP asymmetry obtained from the above distributions 
and projections. 



The technique used to determine the CP asymmetry 
parameters from At distributions is described in detail 
in Ref. [5 . The decay vertex of the signal B meson is 
reconstructed from a kinematic fit of the two D mesons 
to a common vertex including information about the IP 
profile. No information about soft pions is used in the 
vertex reconstruction. The decay vertex and the flavor 
of the accompanying B meson is obtained by an inclu- 
sive approach using the remaining charged tracks that are 
not used in the signal B reconstruction. Requirements 
on the quality of reconstructed B vertices and on the 
number of hits in the silicon vertex detector are applied. 
The algorithms applied to obtain the 6-flavor charge q 
and a tagging quality variable r are described in detail in 
Ref. [35] . The variable r is related to the mistag fractions 
determined from b — > c control samples and ranges from 
r = (no flavor discrimination) to r = 1 (unambigu- 
ous flavor assignment). The data is divided into seven r 
intervals. 

The CP asymmetry parameters are determined by un- 
binned maximum likelihood fits to the At distributions. 
The probability density function used to describe the At 
distributions is given by 

P = (1 - U) Yl fk J P>* W Rk (At - At 1 )} d (At') 



+/01P01 (At) 



(3) 



where the index k denotes signal and background com- 
ponents and the fraction fk depends on the r interval 
and is evaluated on an event-by-event basis as a func- 
tion of Mb c and AE. The signal component consists of 
the convolution of distributions given by modifications of 
Eq. [T] and [2] that include the effect of incorrect flavor as- 
signments and of a resolution function to account for the 
finite resolution of the vertex reconstruction [26 . The 
background component is parameterized by the convo- 
lution of the sum of a prompt and an exponential dis- 
tribution allowing for effective lifetimes and a resolution 
function composed of the sum of two Gaussian functions. 
The parameters of the background components are fixed 
to values determined by fits to Mb c < 5.26 GeV/c 2 side- 
bands. A Gaussian function P Q \ with a broad width of 
about 35 ps and a small fraction f Q \ of about 2 x I0~ 4 is 
added to account for outlier events with large At. 

The free parameters in the B° —> D + D~ fit are S^+jj- 
and Cd+d- an d the free parameters in the B° — > _D* ± D =F 
fit are A D * Dl S D * D , C D * D , AS D * D and AC D * D . The 
lifetime t^o and mass difference Arrid are fixed to current 
world averages [24]. The fits are performed in a signal 
region defined by \AE\ < 30 MeV and 5.27 GeV/c 2 < 
M bc < 5.29 GeV/c 2 . The signal purity is 62% (59%) for 
B° -> D+D- (B° -> D* ± D+). For B° -> D+D- the 
results are 



Sd+d- 
Cd+d- 



-1.06 tofi ±0.08 
-0.43 ±0.16 ±0.05, 



(4) 



G 



TABLE II: Summary of systematic uncertainties in the time-dependent CP asymmetry parameters for B — > D + D and 
B o _^ D *± D T decays (in units of 1CT 2 ). 



Source 


&D+D- 




A 

S\d*D 


OD'D 


^D*D 




A n 

AC D * D 


Vertex reconstruction 


3.6 


2.2 


1.3 


2.5 


2.3 


2.4 


2.3 


At resolution function 


0.0 


2.4 


0.4 


3.0 


l.l 


i n 
1 . 9 


U.O 


Background At PDFs 


2.7 


0.5 


0.2 


0.7 


0.2 


0.5 


0.1 


Signal purity 


1.2 


1.8 


0.2 


0.9 


0.4 


0.3 


0.2 


Physics parameters 


0.7 


0.4 


< 0.1 


0.2 


0.1 


0.2 


< 0.1 


Flavor tagging 


0.7 


0.6 


< 0.1 


0.4 


0.3 


0.3 


0.2 


Possible fit bias 


0.8 


0.2 


0.6 


0.8 


1.1 


0.8 


0.5 


Peaking background 


0.3 


0.9 


0.4 


1.3 


0.5 


0.8 


0.7 


Tag-side interference 


1.4 


3.2 


0.2 


1.1 


3.1 


0.9 


0.6 


Total 


8.2 


5.1 


1.6 


4.9 


4.3 


3.5 


2.6 



and for B° D* ± D zfl 



A D * 


D = 


+0.06 ±0.05 ±0.02 




D = 


-0.78 ±0.15 ±0.05 




D = 


-0.01 ±0.11 ±0.04 


ASd* 


D = 


-0.13 ±0.15 ±0.04 


AC d * 


D = 


+0.12 ±0.11 ±0.03, 



where the first uncertainty is statistical and the second 
systematic. The At distributions and projections of the 
fits are shown in Fig. [2] 

The systematic uncertainties in the CP asymmetry pa- 
rameters are evaluated for each decay mode and are sum- 
marized in Table [IT] Sources of systematic uncertainties 
on the vertex reconstruction are the IP profile constraint, 
requirements on the vertex fit quality for signal and tag- 
ging B mesons, requirements on impact parameters of 
tracks in the reconstruction of the tagging B meson and 
the At fit range. These contributions are estimated by 
variations of each of the applied requirements. Further 
contributions to the vertex reconstruction are a global 
SVD misalignment and a Az bias, which are both esti- 
mated by MC simulations. The contributions due to the 
At resolution functions, the At parameterization of back- 
ground components, the calculation of the signal purity 
and the physics parameters t b o and Am^ are estimated 
by varying the fixed parameters within their uncertain- 
ties. The systematic uncertainty due to flavor tagging is 
estimated by varying the mistag fractions in each r in- 
terval within their uncertainties. A possible fit bias is 
estimated from a large sample of MC simulated signal 
decays. The effect of the peaking background is stud- 
ied using MC simulations allowing for CP violation in 
non-resonant decays. The possible interference between 
Cabibbo-favored b — > cud and suppressed b — > ucd am- 
plitudes in the decay of the tagging B meson, referred to 
as tag-side interference [27], is studied using MC simula- 
tions with inputs obtained from B° — > D*~i + vi control 
samples. The largest deviations in the above MC stud- 
ies are assigned as systematic uncertainties. The total 
systematic uncertainty is obtained by adding all contri- 
butions in quadrature. 



The significance of the results is studied by a 
likelihood-ratio approach. For B° — s- D + D~ we exclude 
the conservation of CP symmetry (Sd+d- = Cd+d- = 
0) at a confidence level of 1 — 2.7 x 10~ 5 corresponding 
to 4.2cr. For B° — > D* ± D T the conservation of CP sym- 
metry (Ad*d = Sd*d = Cd*d = 0) is excluded at a 
confidence level of 1 — 6.8 x 10~ 5 corresponding to 4. Oct. 
These results account for both the statistical and the sys- 
tematic uncertainties. 

The fit procedure was validated by various cross- 
checks. The same analysis was performed for B° — > 
D+DW- decays. The results are A Ds d = -0.01 ± 0.02, 
Sd s d = -0.05 ± 0.05, C d ,d = +0.01 ± 0.03, AS Ds d = 
+0.01 + 0.05 and AC Ds d = -0.95 ±0.03 in B° -» D+D~ 
and A d ,d* = +0.01 ± 0.02, S d ,d* = -0.04 ± 0.05, 
C DsD * = +0.06 ± 0.03, AS DsD * = +0.10 ± 0.05 and 
AC Db d* = -1.00 ± 0.03 in B° -> D+D*~, where the 
uncertainties are statistical only. The results are consis- 
tent with the assumption of no CP violation in B° — s- 
£)+£)(*)- decays. The lifetimes determined by fits to un- 
tagged B° — > D + D~ and B° -» D* ± D T samples are 
consistent with the world average [24] , 

In summary we report measurements of the branch- 
ing fractions and time-dependent CP violating asymme- 
tries in B° -> D + D~ and B° — s- D* ± D T decays using 
the final Belle data sample of (772 ± 11) x 10 6 BB pairs. 
We measure the branching fractions B (_B° — > D + D~) = 
(2.12 ±0.16 ±0.18) x 10~ 4 and B(B°^D* ± D^) = 
(6.14 ±0.29 ±0.50) x 10~ 4 . The measured CP asym- 
metry parameters are S d +d- = —1-06 14 ^ 0-08 and 
C d +d- = -0.43 ± 0.16 ± 0.05 in B° D + D~ and 
A d *d = +0.06±0.05±0.02, S d *d = -0.78±0.15±0.05, 
C d *d = -0.01±0.11±0.04, ASd* d = -0.13±0.15±0.04 
and AC d *d = +0.12 ± 0.11 ± 0.03 in B° D^D*. 
For B° — > D + D~, the CP asymmetries are approx- 
imately 0.5cr outside of the physical parameter space 

defined by \Js d+d _ +C^ +£) _ < 1 and the direct CP 
asymmetry deviates from zero by approximately 2. Oct. 
For B° D^DT, if the contribution of penguin di- 
agrams is negligible and if the hadronic phase between 
B" ->■ D*+D- and B° ->■ D*~ D+ amplitudes is zero 
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and their magnitudes are the same, then Ad*d, Cd*d, 
A(Sd«d and ACd*d vanish and Sd*d is equal to sin20i. 
Our result is consistent with the above and we measure 
sin 20i = -0.78 ± 0.15 ± 0.05. The CP asymmetries 
obtained in B° — > D + D~ and B° — > D* ± D T decays are 
both in agreement with measurements of decays involving 
b — > (cc)s transitions [H [S] and with previous measure- 
ments of B° ->■ L»(*)±d(*)=f decays 0HDHI2]. We find 
evidence for CP violation in both decay channels with 
a significance of >4er. These results supersede previous 
measurements of branching fractions and time-dependent 
CP asymmetries in B° -> D+D~ and B° -> D* ± D =F by 
the Belle Collaboration [51 ITUl 135] . 
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